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ABSTRACT 

We report the detection of CN(7V=3^2) emission towards the Cloverleaf quasar {z = 2.56) based 
on observations with the IRAM Plateau de Bure Interferometer. This is the first clear detection of 
emission from this radical at high redshift. CN emission is a tracer of dense molecular hydrogen 
gas (n(H2) > lO^cm"^) within star-forming molecular clouds, in particular in regions where the 
clouds are affected by UV radiation. The HCN/CN intensity ratio can be used as a diagnostic for 
the relative importance of photodissociation regions (PDRs) in a source, and as a sensitive probe 
of optical depth, the radiation field, and photochemical processes. We derive a lensing-corrected 
CN(Ar=3^2) line luminosity of i^cN(3-2) = (^-^ ^ ^•^) ^ ^ kms^^pc^ The ratio between CN 
luminosity and far-infrared luminosity falls within the scatter of the same relationship found for \ow-z 
(ultra-) luminous infrared galaxies. Combining our new results with CO(J=3^2) and HCN(J=1^0) 
measurements from the literature and assuming thermal excitation for all transitions, we find a CO/CN 
luminosity ratio of 9.3 ± 1.9 and a HCN/CN luminosity ratio of 0.95 =b 0.15. However, we find that 
the CN(A^=3^2) line is likely only subthermally excited, implying that those ratios may only provide 
upper limits for the intrinsic 1^0 line luminosity ratios. We conclude that, in combination with other 
molecular gas tracers like CO, HCN, and HCO+, CN is an important probe of the physical conditions 
and chemical composition of dense molecular environments at high redshift. 

Subject headings: galaxies: active, starburst, formation, high redshift — cosmology: observations — 
radio lines: galaxies 
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1. INTRODUCTION 

Investigations of the dense molecular interstellar 
medium (ISM) in high-redshift galaxies are of fundamen- 
tal importance to further our understanding of the early 
phases of galaxy formation and evolution, as it harbors 
the environments in which the actual star formation is 
believed to occur. Recent studies of the molecular gas 
phase in high-z galaxies using CO emission lines have 
revealed large molecular gas reservoirs with masses in ex- 
cess of 10^^ M0 (see review by Solomon & Vanden Bout 
I2005L and references therein). 

However, while the relative brightness of CO emission 
lines renders this molecule the most common tracer of 
the requisite material for star formation, the low den- 
sities of only ^ 10^ — 10^ cm~^ required to excite 
its lower-J transitions (due to its low dipole moment of 
only = 0.11) imply that CO is not a specific tracer 
of dense molecular cloud cores, i.e., the regions where 
stars are actively formed. 

In contrast, recent studies of nearby actively star- 
forming galaxies have shown that the high dipole 
molecule HCN (/ig^^ = 2.98) is a far better tr acer of 
such dense molecular cores (e.g. Gao & Solomon l2004al . 
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I2004b[ ). The critical density of nu^ ^ 10^ - 10^ cm-^ 
to collisionally thermalize its lower-J transitions is much 
higher than that of CO, and of the same order as the den- 
sities found in bright star-forming regions in the Galaxy 
(e.g. the Orion Bar, Hogerheijde et al. 1995). A main re- 
sult in this context is the finding that the HCN luminos- 
ity correlates well with the far-infrared (FIR) luminosi- 
ties (which is commonly used to estimate star-formation 
rates at high redshift) over 7-8 orders of magnitude, from 
Galactic dense cores to the highest redshift quasars (Wu 
et al. 2005). 

However, in nearby luminous and ultra-luminous in- 
frared galaxies (LIRGs/ULIRGs), it has been found that 
systems with similar HCN/CO and HCN/FIR luminos- 
ity ratios may have quite different dense gas properties 
regarding their chemical composition, and their gas ex- 
citation (e.g., Aalto et al. 2002). To better understand 
the physical and chemical state of the dense molecular 
gas phase which directly relates to star formation, it 
has proven essential to study additional bright tracers 
of dense gas with properties different from HCN, such 
as the cyanide radical (CN). Due to its lower dipole mo- 
ment relative to HCN (/i§^ = 1.45), its critical den- 
sity is lower by about a factor of 5. Observations of 
CN emission toward the Orion A molecular cloud com- 
plex have shown that CN filaments trace the dense inter- 
faces between the molecular cloud and the major ion- 
ization fronts (Rodriguez-Franco et al. 11998). It has 
also been found that the [CN]/[HCN] abundance ratio 
is greatly enhanced in the central region of the star- 
burst galaxy M82, beingas high as ~5 across the entire 
nucleus (Fuente et al. l2005f ). These observations indi- 
cate that CN is a good tracer of gas layers which are 
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affected by photochemistry, since this molecule appears 
to be predominantly found in regions exposed to ioniz- 
ing stellar UV radiation. It has also been found that, 
due to the rapid destruction of other dense gas tracers 
like HCN, the abundance of CN tends to be enhanced in 
areas where the UV radiation field is only partly atten- 
uated, such as in zones close to the sur face o f photodis- 
sociation regions (PDRs, Fuente et al. 1 19931 : Sternberg 
& Dalgarno 1995; Jansen et al. 1995). The CN/HCN 
intensity ratio can thus be used as a diagnostic for the 
relative importance of PDRs in a source, and a sensitive 
probe of optical depth, the radiation field, and p hoto- 
chemical processes (e.g., Boger & Sternberg '2005"). CN 
emission may thus in a sense be a more specific tracer 
for star formation than, e.g., HCN, which only traces re- 
gions of dense gas in general. In addition, theoretical 
studies of the chemical composition of molecular gas for 
different ionization parameters suggest that the relative 
abundance of CN may also be enhanced in X-ray domi- 
nated regions (XDRs), such as AGN environments (Lepp 
& Dalgarno 1996). This exemplifies why it is desirable to 
search for CN emission, which is part of our current ef- 
fort to study molecular tracers other than CO out to high 
redshifts. Such an investigation is imperative to obtain 
more meaningful constraints on the physical properties 
and chemical composition of the dense molecular ISM in 
distant galaxies. 

In this paper, we report the first high-z detection of 
CN(A/'=3^2) emission, which was observed towards the 
Clover leaf quasar {z = 2.56) with the IRAM Plateau de 
Bure Interferometer (PdBI)^. Due to its strong gravita- 
tional magnification (magnification factor /iL = H, Ven- 
turini & Solomon ^03), the Cloverleaf is the brightest 
CO source at high redshift (e.g., Barvainis et al. 1994), 
and one of the most prolific sources of molecular lines 
beyond z=2. It was the first z > 2 source to be de- 
tected in HCN (Solomon et al. 2003) and HCO+ (Riech- 
ers et al. ■2006 a) emission. It was also detec ted in both 
C I fine stru cture lines (Barvainis et al. Il997l : WeiB et al. 
I2003L l2005l ) . We use a standard concordance cosmology 
throughout, with Hq = 71 kms~^Mpc~\ flu = 0.27, 
and Qa = 0.73 (Spergel et al. 120031 . 12007[ ). 

2. OBSERVATIONS 

We observed the CN(A^=3^2) transition line (i^rest = 
339.4467770 - 340.2791661 GHz for the different fine 
structure [fs] and hyperfine structure [hfs] transitions 
of the v=0 vibrational state) towards H1413+117 (the 
Cloverleaf quasar) using the PdBI in D configuration 
between 2006 July 26 and September 03. A t the tar- 
get redshift of z=2. 55784 (WeiB etal. 120031 ). the line 
is shifted to ~95.6GHz (3.14mm). The total integra- 
tion time amounts to 17.5 hr using 5 antennas, result- 
ing in 5.8 hr equivalent on-source time with 6 antennas 
after discarding unusable visibility data. The nearby 
sources 1354+195 and 1502+106 (distance to the Clover- 
leaf: 9.0° and 12.0°) were observed every 20 minutes for 
pointing, secondary amplitude and phase calibrations. 
For primary fiux calibration, several nearby calibrators 
(MWC349, CRL618, 3C273, 3C345, and NRAO 150) 
were observed during all runs. 

6 IRAM is supported by INSU/CNRS (France), MPG (Ger- 
many), and IGN (Spain). 



Vobs [GHz] 
95.8 95.7 95.6 95.5 95.4 



0.5 




Model spectrum of CN(N=3-2) in the Cloverleaf 


/ 3 












- 7 









340.8 340.4 340 339.6 339.2 
Vrest [GHz] 



Fig. 1. — Model spectrum of the CN(A^=3^2) emission in the 
Cloverleaf. The horizontal axes indicate the rest and observed fre- 
quencies. The full range of the PdBI bandpass is shown. The 
vertical axis indicates the predicted intensity, normalized to the 
(Ar=3^2, J=|^|, F=|^|) component. The relative intensi- 
ties are computed for the LTE case. The crosses indicate the hy- 
perfine structure components. The plus signs indicate the summed 
intensities of the components in three different frequency bins. The 
dashed lines are Gaussian fits to the summed intensities, assuming 
the width of the CO( J=3^2) fine (WeiB et al. 2003) for the sub- 
components (labeled 'l'-'3'). The solid line is a sum of all three 
components, and indicates the model-predicted CN(A^=3^2) line 
shape. The thick dotted line is a single Gaussian (used to fit the 
observations), fitted to the model line. The dotted vertical line 
indicates the tuning frequency, corresponding to zero velocity in 
Fig.H] 

The correlator was tuned to a frequency of 95.603 GHz, 
which corresponds to the central position between the 
brightest hfs transitions^ of the CN(A/'=3^2) line at 
95.6332 and 95.5736 GHz. The total bandwidth of 
580 MHz (~1800 kms~^) used for the observations is 
large enough to cover all fs and hfs transitions of the 
CN(A^=3^2) line. It also covers enough channels that 
are free of line emission to constrain the 3 mm continuum 
emission of the Cloverleaf. 

For data reduction and analysis, the IRAM GILD AS 
package was used. All data were mapped using the 
CLEAN algorithm and 'natural' weighting without ap- 
plying a further taper; this results in a synthesized beam 
of 5.3^^x4.9 (-42kpc at z = 2.56). The final rms in the 
combined map is 0.25 mJy beam~^ for a 190 MHz (cor- 
responding to 596 kms~^) channel, 0.4 mJy beam" for 
a 80 MHz (252 kms~^) channel, and 0.8 mJy beam~^ for 
a 20 MHz (63 kms"^) channel. 

3. THEORETICAL CONSIDERATIONS 

Due to its fine stucture and hyperfine structure split- 
ting, the CN(A/'=3^2) emission line is distributed over 
19 lines^, in 3 main components separated by more than 

^ These components are actually blends of the [A/'=3^2, 
J=|^|, l^f, and f^|)] hfs components at 

i/rest = 340.0315440 - 340.0354080 GHz and [A^=3^2, J=|^|, 

, and f^f)] hfs components at z^rest 
W 10.2 177700 - 340.2485764 GHz (components 2 and 3 in Fig.[l]). 

^ We assume that CN is in its ground electronic state (^XI), and 
that the spins couple according to Hund's case (b) coupling scheme: 
N-\-S = J, (fs coupling) and J + /c+/n = ^ (hfs coupling). Here, 
N is the rotational angular momentum vector, S is the electronic 
spin, and / is a nuclear spin. 
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200 MHz from each other in the rest frame. This separa- 
tion is of the same order as the kinematical broa denin g; 
of the CO hnes in the Cloverleaf (WeiB et al. 120031 ). 
causing the hfs components to be blended. To analyze 
the intrinsic line shape of the CN(A^=3^2) transition in 
the Cloverleaf, we calculated a synthetic line profile, as- 
suming optically thin emission in Local Thermodynamic 
Equilibrium (LTE) to derive the relative intensities of 
the hfs components. The relative intensities of the com- 
ponents were computed using the laboratory data from 
Skatrud et al. (1983), and app roximation (4) of equa- 
tion (1) of Pickett et al. (|l998f ). Assuming that the 3 
main CN(A/'=3^2) components are kinematically broad- 
ened in the same way as the C O lines (416 ± 6 kms~^ 
FWHM; see WeiB et al. l2003f ). we obtain the synthetic 
line profile displayed in Fig. [T] (solid line). Under the 
given assumptions, the contribution from all hfs lines of 
component 1 to the total intensity are negligible, and 
not detectable at the given signal-to-noise (see below). 
A single Gaussian fit (dotted line, 484 kms~^ FWHM, 
or 116% of the CO lines) agrees with the more detailed 
model profile (solid line) within a few percent in its peak 
position and integral, i.e., well within the observational 
errors. Note however that the width and peak position 
of this Gaussian depends on the relative intensities of 
the hfs components, which in turn depend on the above 
assumptions, in particular the optical depth. In the op- 
tically thin LTE case, components 2 and 3 have a peak 
strength ratio of 1:1.4. As an example, if they had a ra- 
tio of unity, the fitted Gaussian to the line profile would 
have a FWHM of 603 kms-\ or 145% of the CO lines. 
We thus conclude that fitting a Gaussian to the observed 
line profile contains all relevant information, while min- 
imizing the number of free fit parameters, and thus is 
preferred over a more complex fitting procedure to de- 
scribe and analyze the CN(A^=3^2) profile at the given 
signal-to-noise ratio. 

4. RESULTS 
4.1. CN Maps and Spectrum 

We have detected emission from the CN(A/'=3^2) 
transition line towards the Cloverleaf quasar {z = 2.56). 
The velocity-integrated CN(A^=3^2) line map is shown 
in Fig. O The cross indicates the geometrical cen- 
ter position of the resolved CO(J=7^6) map at a = 
14hi5m4gs 233^ (5 = +11°29'43''.50 (Alloin etal. 1997). 
The line emission is clearly detected at 7cr over a range 
of 596 kms~^ (190 MHz); the source appears unresolved. 
In Fig.[3l six channel maps (252 kms~^, or 80 MHz each) 
of the CN(A^=3^2) line emission are shown. At an rms 
of 0.4mJy beam~^, the line is detected at 4 and 5cr in 
the central channels, and the decline of the line intensity 
towards the line wings is clearly visible toward the outer 
channels, as expected. 

Fig. [4] shows the spectrum of the CN(A/'=3^2) emis- 
sion at a resolution of 63 kms~^ (20 MHz). Zero velocity 
corresponds to the tuning frequency of 95.603 GHz. The 
solid line shows a Gaussian fit to the spectrum. From the 
fit, we derive a line peak flux density of 1.94 ± 0.24 mJy 
beam-\ and a line FWHM of 666 ± 97 kms-^ The fit 
provides an upper limit for the continuum emission at 
the line frequency of <0.25mJy. From the channels as- 
sumed to be free of line emission, we derive a formal 3<j 
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Fig. 2. Velocity-integrated PdBI map of CN(A^=3^2) emis- 
sion towards the Cloverleaf quasar. At a resolution of 5.3^^x4.9^' 
(as indicated in the bottom left corner), the source is unresolved. 
The cross indicates the geometrical center of the CO emission in 
the Cloverleaf (Alloin et al. 1997, see text). Contours are shown at 
(-3, -2, 2, 3, 4, 5, 6, 7)xa (la = 0.25 mJy beam"!). 

upper limit to the continuum peak flux density of 0.7 mJy 
beam~^. Note that the model fit to the dust SED of the 
Cloverleaf by WeiB et al. (2003) would suggest a contin- 
uum flux of ~0.3 mJy at the CN(A^=3^2) line frequency, 
which agrees within the errors with the above estimates. 
As the continuum however is not detected, we do not 
subtract a continuum component from the observed spec- 
trum. This leads to an integrated CN(A/'=3^2) line flux 
of 1.37 ± 0.17Jy kms~^. The velocity offset of the 
Gaussian peak relative to the tuning frequency is -105 
± 38 kms-^ 

The FWHM of the Gaussian fit to the CN(7V=3^2) 
emission line suggests that its components 2 and 3 have 
similar peak strengths (different from 1:1.4 as predicted 
for the optically thin LTE case; Fig. [T]), which may in- 
dicate that the emission is optically thick. Within the 
limited signal-to-noise of the spectrum, the shape of the 
emission line remains compatible with a two-component 
structure. It may thus be possible that the two brightest 
hfs complexes (components 2 and 3) are detected indi- 
vidually, but observations at higher signal-to-noise are 
needed to confirm this result. In any case, a double Gaus- 
sian fit to the line profile gives the same integrated line 
flux as the single Gaussian fit within the errors. 

4.2. Line Luminosities 

From our observations, we derive a CN(A/'=3^2) line 
luminosity of ^cn(3_2) = (4.5 ± 0.5) x 10^ K km s~^ pc^ 
(corrected for gravitational magnification, /iL = H; see 
Table 1 for details). This corresponds to 11% of ^co(3-2) 
(WeiB et al. I2003D . 105% of i^HCN(i-o) (Solomon et 
al. I2003D . and 129% of i^HCO+(i-o) (Riechers et al. 
l2006a[ ) . Assuming that CN is thermally excited up to the 
A/'=3^2 transition (i.e., ^cn(3_2) = ^cn(i-o)) ^^^^ 
optical depth effects can be neglected, this would mean 
that CN emission in the Cloverleaf is slightly brighter 
than that from the other high density probes HCN and 
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TABLE 1 

Molecular line luminosities in the Cloverleaf. 







[mJy] 


[lO^K kms-ipc^] 


Ref. 


CN(Ar=3- 


-2) 


1.94 ± 0.24 


4.5 ± 0.5 


1 


HCN(J=1- 


-0) 


0.24 ± 0.04 


4.3 ± 0.5 


2 


HCO+(J=l- 


-0) 


0.19 ± 0.03 


3.5 ± 0.3 


3 


CO(J=3- 


-2) 


30 ± 1.7 


42 ± 7 


4 



Fig. 3. — Channel maps of the CN(A/'=3^2) emission (same region is shown as in Fig. [2]). One channel width is 80 MHz, or 252 km s ^ 
(at 95.4355, 95.5155, 95.5955, 95.6755, 95.7555, and 95.8355 GHz; frequencies increase with channel number). Contours are shown at (-3, 
-2, 2, 3, 4, 5)xcr (Icr = 0.4 mJy beam" -'^). The beam size (5.3^^x4.9^^) is shown in the bottom left corner; the cross indicates the same 
position as in Fig. [2] 

HCO+. This is remarkable, as HCN and HCO+ have 
higher critical densities than CN in the ground-state 
transition, but the critical density of CN(7V=3^2) is 
higher than that of HCN(J=1^0) and HCO+(J=1^0). 
This may indicate a higher relative filling factor of CN 
relative to HCN and HCO+, or even a relatively high 
chemical abundance of CN. However, with the observa- 
tions existing at present, it is not possible to disentangle 
excitation effects from chemical effects. 

In addition, observations of HCN(J=4^3) to- 
wards the Cloverleaf quasar have shown that 
^HCN(4-3)/^HCN(i-o) — ^'^^ (Solomou et al. l2QQ3l . cap- 
tion of their Tab. 1, and M. Guelin 2007, priv. comm.). 
This ratio is significantly lower than 1, which implies 
that the 4^3 transition of HCN is sub-thermally 
excited. As the CN(A^=3^2) transition has a some- 
what lower but comparable critical density relative to 
HCN(J=4^3), it is likely that the CN(7V=3^2) tran- 
sition is also sub-thermally excited, and the 'intrinsic', 
thermalization-corrected ^cn(i-o) ^^^^ even higher. 
This would be consistent the finding that CN is clearly 
subthermally excited in nearby LIRGs and ULIRGs 
(Aalto et al. l2002f ). This would also imply that CN is 
a brighter tracer of dense gas than HCN and HCO+ in 



References. — [1] This work, [2] Solomon etal. (120031) . [3] Riechers 
et al. (2006a), [4] Weifi etal. (2003). 

Note. — Luminosities are corrected for gravitational magnification. 

this high redshift quasar. 

5. DISCUSSION 

In the following, we discuss various relationships be- 
tween the emission observed in CN and other molecules 
(CO and HCN) and the far-IR continuum for a sam- 
ple o f local spiral/starburst galaxies (Henkel et al. Il988l . 
119981 . 2007, priv. comm.; Wang et al. 2004; additional 
CO/HCN/ FIR d ata from Eckart et al. 1990; Nguyen-Q- 
Rieu et al. 1 19921 : Mauersberger et al. 1996; Mao et al. 
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Fig. 4. — Spectrum of the C'N(N=3—^2) emission at a resolution 
of 63 kms"-*^ (20 MHz). The velocity scale is relative to the tuning 
frequency of 95.603 GHz. The rms per velocity bin is 0.8 mJy. The 
solid line shows a Gaussian fit to the data. 



2000; Sanders et al. 2003), low-z (U)LIRGs (Aalto et al. 
2002), and the C loverleaf (this work; Solomon et al. 2QQ3I : 
WeiB et al. l2003D as shown in Fig. El Note that, due to the 
weak signal-to-noise ratio of Icr, we do not include the 
CN(7V=4^3) observations of APM 08279+5255 (Guelin 
et al. 2003, their Tab. 1) in this discussion. 

As CO (J= 1—^0) is not yet detected toward the Clover- 
leaf, we assume in the following that CO is thermalized 
up to the 3^2 transition, i.e., ^co(i-o) ~ ^co(3-2) (^^^ 
Riechers et al. l 2006bl for justification). Both effects of 
thermal and sub-thermal excitation of the CN(A^=3^2) 
emission line in the Cloverleaf are discussed. CO, CN, 
and HCN line intensities for all low-2: galaxies were 
measured directly in the ground-state transitions, all 
quoted luminosities thus are 1^0 luminosities. We do 
not discuss effects of differential lensing^, as models in- 
dicate similar sizes for molecular and dust emission in 
the Cloverleaf (Solomon et al. [2003). 

Figure 5a: Wu et al. (|2005l ) find a linear correlation 
between i^HCN ^fir over 7-8 orders of magnitude in 
luminosity, which is believed to indicate that the dense 
gas tracer HCN is also a good estimator for star forma- 
tion rates out to high z. If CN were to be a good tracer 



^ Note that Chartas et al. (|2007l ') recently reported a likely mi- 
crolensing event toward the Cloverleaf; however, this event took 
place close to 2000 April, so 3-7 years before all of the molecular 
gas observations of this source discussed in this paper (obtained 
between 2003 March and 2006 September). 



CN(7V=3-2) in the Cloverleaf 



alO^ 







109 



: 






Mrk231 1— @H 


_ 


Arp220 IQj/ 




/ie|NGC6240 - 






: M82 / 




As 

IC342 / NGC253 

- 




1 I'll 


^spirals/ starbursts- 
OLIRGs/ULIRGs 
1 X^loverleaf QSO 



blO^ 



10^ 



© 



106 10^ 108 109 

L'cN pc^] 



IQIO 



109 



H— K 



NGC6240 O 
O 



NGC253 -jP 



^ spirals/starbursts 
O LIRGs/ULIRGs 
^ XCIpverleaf QSO ^ 



ClO^ 



10^ 





^10^ 



10 



20 
d/^ cn 



30 



109 



z 1 1 1 1 1 1 


Mrk231 
® 


A 220O 




- o ® 

o 


NGC6240O - 




o - 


: M82 




NGC253W 




< < , , 1 < 


IC342 

i:^ spirals/starbursts- 
O LIRGs/ULIRGs 
XCIoverleaf QSO 



Fig. 5. — CN luminosity relations for a sample of local spiral/starburst galaxies, low-z IR-luminous galaxies, and the Cloverleaf. The 
Cloverleaf luminosities are corrected for gravitational lensing (/iL = 11)- The two symbols with star insets denote tentative CN detections. 
The solid line is a least squares fit to all data except the Cloverleaf. The error bars indicate the statistical errors of the line luminosity 
measurements. The excitation of the CN(A^=3^2) transition in the Cloverleaf may be sub-thermal, which would affect the (extrapolated) 
CN(A^=1^0) line luminosity as indicated by the horizontal, shaded regions. As an example, the vertical bars on the shaded regions show 



a case where 



CN(3-2) 



: 0.5 X 



CN(l-O) 



in the Cloverleaf. See text for more details. 



of dense, actively star-forming g as, a similar trend may 
be expected between L'^^ and Lfir- Figure [5^1 shows 
that correlates closely with Lpm; a linear least 

squares fit to all galaxies (excluding the Cloverleaf) yields 
log(LFiR) = (0.89 ± 0.09) X log(L^N) + (4-2 ± 0.7). This 
slope may suggest a trend of decreasing CN/FIR lumi- 
nosity ratio towards higher luminosities (see also Aalto 
et al. I2OO2I ). It however also is consistent with unity 
within the statistical uncertainties (which are quite large 
due to the limited sample size) and systematical errors 
(e.g., different beam sizes for CN line and IR continuum 
observations), and thus with CN being a valuable tracer 
of star formation. The Cloverleaf agrees remarkably well 
with the relation defined by the \ow-z galaxies (even if 
the CN(A^=3^2) transition is sub-thermally excited, as 
indicated by the horizontal bar on the shaded region), 
extending the observed trend to higher luminosities and 
out to high redshift. 

Figure 5b: Gao & Solomon (|2004bl ) observe a rise in 
dense gas fraction with LpiR. The Lqq/Lq^ ratio may 
be considered (the inverse of) a tracer of the dense gas 
fraction. A significant spread in the Lqq/Lq^ ratio is 
found for the sample shown in Fig. [6)3, ranging from a 
few to almost a factor of 30. There is no obvious trend 
of Lqq/Lq^ with I/fir; however, the sample size is too 
small and heterogenous to come to a definite conclusion. 
Note that, due to the likely sub-thermal excitation of 
the CN(A^=3^2) transition (as indicated by the shaded 
region), the Lqq/L^^ ratio may be considered an upper 
limit of the 'intrinsic', thermalization-corrected ratio of 
the 1^0 transitions for the Cloverleaf. 

Figure 5c: Assuming that CN traces the UV field 
(which falls off rapidly with growing distance from star- 
forming regions, or an active galactic nucleus), while 
HCN traces dense gas in general, it would be expected 
that the HCN/CN ratio stays constant or decreases with 
increasing LpiR. In this scenario, a decrease in the 
HCN/CN ratio would correspond to an increase in the 
filling factor of UV-illuminated clouds with Lfir. This 
would however also cause a rise in L^js^/LpiR, which 
is not observed (Fig. [5^). If differences in the rela- 
tive chemical abundances of CN and HCN and optical 



depth/thermalization effects do not play a major role, 
one may thus expect ^hcn/^cn ^^^y ^^^^^Y constant 
with LpiR. Based on the observations of their subsample, 
Aalto et al. (2002) suggest that the HCN/CN intensity 
ratio may increase slightly with LpiR. Such a trend how- 
ever is not seen in the larger sample displayed in Fig. [5t 
(note that ^hcn/^cn Cloverleaf may again be 

considered an upper limit due to possible sub-thermal 
excitation of CN(A^=3^2), as indicated by the shaded 
region). Clearly, improved statistics and better models 
are required to address this issue in more detail. 

Chemical models suggest that CN is produced quite ef- 
ficiently in the presence of a strong UV field, which also 
leads to an enhanced ionization rate (Boger & Stern- 
berg 2^. The large observed [CN]/[HCN] abundance 
ratio of ^5 toward the molecular disk of M82 may be 
indicative of a large dense PDR bathed in the intense 
radiation field of the starburst environment (Fuente et 
al. 2005), lending observational support to these mod- 
els. Other models suggest that the relative abundance 
of CN may also be enhanced in regions with elevated X- 
ray ionizati on rat es, such as AGN environments (Lepp 
& Dalgarno 119961 : Meijerink & Spaans 2005). Especially 
in the Cloverleaf, where the central quasar is known to 
contribute a significant fraction to the heatin g of t he gas 
and dust (WeiB et al. 2003; Solomon et al. l2003f ). such 
a scenario would appear reasonable and may explain the 
relatively high brightness of CN compared to other dense 
gas tracers in the Cloverleaf. However, note that in 
Fig. 5c, the pure starburst galaxies show the brightest 
CN emission relative to HCN, and galaxies with a rela- 
tively strong AGN like Mrk 231 and NGC 6240 have some 
of the highest ^hcn/^cn i"atios. This may indicate that 
the strength of the UV field in a starburst environment 
has a significantly greater impact on the global gas-phase 
production rate of CN in a galaxy than the strength of 
the X-ray field, lending further support to the supposi- 
tion that CN is a good star formation tracer. We thus 
conclude that the XDR scenario alone does not explain 
the observed dense gas properties of the (U)LIRGs shown 
in Fig. 5 without any further assumptions. 

Overall, the Cloverleaf follows the L^js^-Lfir relation 
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as defined by nearby galaxies of different types remark- 
ably well, and extends this relation to higher luminosi- 
ties. This relation now appears to hold over almost 3 
orders of magnitude. If these findings were to hold for 
other distant quasars, CN would be an excellent probe 
of dense molecular environments out to high redshifts. 
Also, the higher-order transitions of CN appear to be 
brighter than those of other dense gas tracers like HCN. 
This may prove to be of particular importance for fu- 
ture high-z studies with the Atacama Large Millime- 
ter /submillimeter Array (ALMA), which will offer the 
opportunity to probe to fainter galaxy populations in 
general, but will also be restricted to the higher-order 
transitions of the most common dense molecular gas trac- 
ers (typically 3^2 and higher at z > 2). 



with the IRAM Plateau de Bure Interferometer. IRAM 
is supported by INSU/CNRS (France), MPG (Ger- 
many), and IGN (Spain). D. R. acknowledges sup- 
port from the Deutsche Forschungsgemeinschaft (DFG) 
Priority Programme 1177. C. C. acknowledges sup- 
port from the Max-Planck- Gesellschaft and the Alexan- 
der von Humboldt- Stiftung through the Max-Planck- 
Forschungspreis 2005. The authors would like to thank 
Susanne Aalto for fruitful discussions on the subject mat- 
ter. We also would like to thank the staff at IRAM, in 
particular Jan Martin Winters, for their assistance in set- 
ting up and carrying out the observations, as well as for 
providing a preliminary reduction of the data. The au- 
thors would like to thank the referee for useful comments 
that helped to improve the manuscript. 



This research is based on observations carried out 



REFERENCES 



Aalto, S., Polatidis, A. G., Hiittemeister, S., & Curran, S. J. 2002, 
A&A, 381, 783 

Alloin, D., Guilloteau, S., Barvainis, R., Antonucci, R., & Tacconi, 

L., 1997, A&A, 321, 24 
Barvainis, R., Tacconi, L., Antonucci, R., Alloin, D., &: Coleman, 

P. 1994, Nature, 371, 586 
Barvainis, R., Maloney, P., Antonucci, R., &: Alloin, D. 1997, ApJ, 

484, 695 

Boger, G. I., & Sternberg, A. 2005, ApJ, 632, 302 

Chartas, G., Eracleous, M., Dai, X., Agol, E., &; Gallagher, S. 2007, 

ApJ, in press (astro-ph/0702742 ) 
Eckart, A., Downes, D., Genzel, R., Harris, A. I., Jaffe, D. T., & 

Wild, W. 1990, ApJ, 348, 434 
Fuente, A., Mart in- Pint ado, J., Cernicharo, J., &: Bachiller, R. 

1993, A&A, 276, 473 
Fuente, A., Garcia-Burillo, S., Gerin, M., Teyssier, D., Usero, A., 

Rizzo, J. R., & de Vicente, P. 2005, ApJ, 619, L155 
Gao, Y., & Solomon, P. M. 2004a, ApJS, 152, 63 
Gao, Y., & Solomon, P. M. 2004b, ApJ, 606, 271 
Guelin, M., Salome, P., Neri, R., Garcia-Burillo, S., Gracia-Carpio, 

J., et al. 2007, A&A, 462, L45 
Henkel, C, Schilke, P., & Mauersberger, R. 1988, A&A, 201, L23 
Henkel, C, Chin, Y.-N., Mauersberger, R., &; Whiteoak, J. B. 1998, 

A&A, 329, 443 

Hogerheijde, M. R., Jansen, D. J., Sz van Dishoeck, E. F. 1995, 
A&A, 294, 792 

Jansen, D. J., van Dishoeck, E. F., Black, J. H., Spaans, M., &: 

Sosin, C. 1995, A&A, 302, 223 
Lepp, S., &; Dalgarno, A. 1996, A&A, 306, L21 
Mao, R. Q., Henkel, C, Schulz, A., Zielinsky, M., Mauersberger, 

R., Storzer, H., Wilson, T. L., Gensheimer, P. 2000, A&A, 358, 

433 

Mauersberger, R., Henkel, C, Wielebinski, R., Wikhnd, T., &; 
Renter, H.-P. 1996, A&A, 305, 421 



Meijerink, R., &; Spaans, M. 2005, A&A, 436, 397 
Nguyen-Q-Rieu, Jackson, J. M., Henkel, C, Truong-Bach, & 

Mauersberger, R. 1992, ApJ, 399, 521 
Pickett, H. M., Poynter, R. L., Cohen, E. A., Delitsky, 

M. L., Pearson, J. C, &; Miiller, H. S. P. 1998, 

J. Quant. Spectrosc. Radiat. Transfer, 60, 883 
Riechers, D. A., Walter, F., CariUi, C. L., WeiB A., Bertoldi, F., 

Menten, K. M., Knudsen, K. K., & Cox, P. 2006a, ApJ, 645, L13 
Riechers, D. A., Walter, F., CariUi, C. L., et al. 2006b, ApJ, 650, 

604 

Rodriguez-Franco, A., Martin-Pintado, J., &; Fuente, A. 1998, 
A&A, 329, 1097 

Sanders, D. B., Mazzarella, J. M., Kim, D.-C, Surace, J. A., & 

Soifer, B. T. 2003, AJ, 126, 1607 
Skatrud, D. D., de Lucia, F. C, Blake, G. A., & Sastry, K. V. L. N. 

1983, J. Mol. Spectrosc, 99, 35 
Solomon, P., Vanden Bout, P., CariUi, C, &; GueUn, M. 2003, 

Nature 426, 636 
Solomon, P. M., & Vanden Bout, P. A. 2005, ARA&A, 43, 677 
Spergel, D. N., Verde, L., Peiris, H. V., et al. 2003, ApJS, 148, 175 
Spergel, D. N., Bean, R., Dore, O., et al. 2007, ApJ, in press 

(astro-ph/0603449) 
Sternberg, A., &; Dalgarno, A. 1995, ApJS, 99, 565 
Venturini, S., &; Solomon, P. M. 2003, ApJ, 590, 740 
Wang, M., Henkel, C, Chin, Y.-N., et al. 2004, A&A, 422, 883 
WeiB, A., Henkel, C, Downes, D., &; Walter, F. 2003, A&A, 409, 

L41 

WeiB, A., Downes, D., Henkel, C, Sz Walter, F. 2005, A&A, 429, 
L25 

Wu, J., Evans H, N. J., Gao, Y., Solomon, P. M., Shirley, Y. L., & 
Vanden Bout, P. A. 2005, ApJ, 635, L173 



